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We discuss the finite-temperature phase diagram in the three-dimensional Bose-Hubbard (BH) 
model in the strong correlation regime, relevant for Bose-Einstein condensates in optical lattices, 
by employing a quantum rotor approach. In systems with strong on site repulsive interactions, the 
rotor U(l) phase variable dual to the local boson density emerges as an important collective field. 
After establishing the connection between the rotor construction and the the on-site interaction 
in the BH model the robust effective action formalism is developed which allows us to study the 
superfluid phase transition in various temperature-interaction regimes. 
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I. INTRODUCTION 



The ability to confine ultracold quantum gases in op- 
tical lattices is already having a major impact in fields 
as diverse as condensed-matter physics and quantum in- 
formation processing's An optical lattice is essentially 
a periodic intensity pattern that is formed by the inter- 
ference of two or more laser beams. The simplest opti- 
cal lattice consists of the region that is formed when two 
laser beams with the same wavelength travelling in oppo- 
site directions meet each other and form an interference 
pattern. This artificial structure is able to trap an atom 
because the electric fields of the lasers induce an elec- 
tric dipole moment in the atom. However, the quantum- 
mechanical tunnelling allows the atoms to spread through 
the optical lattice to some degree. In a Bose-Einstein con- 
densate, this tunnelling process dominates the behaviour 
of the atoms, which causes the system to have a per- 
fect phase coherence between the matter waves on differ- 
ent lattices sites. This quantum phase transitions of the 
Bose-Einstein condensates loaded into the lowest vibra- 
tional level of single wells of an optical lattice in the strict 
sense can exist only at temperature T = 0. However, in 
typical experimental situations we must take into consid- 
eration thermal fluctuations in the particle number per 
site. In the presence of the finite-temperature the nonzero 
value of the compressibility is expected in contrast to 
the incompressible Mott state at T = 0. The experimen- 
tal data only signals that the system nears a quantum 
phase transition if the temperature is extrapolated to 
zero. What the experiments really observe is a transition 
from the superfluid to the normal liquid whose compress- 
ibility is very close to zero and the system is practically 
a Mott insulator. This issue has gained recently much 
attentiorPHll. TheoreticaP and numericaPm' approaches 
have designed to these systems in three-dimensions but 
only recently the finite-temperature effects have been 
studied systematical!}^^. Yu at. alP^ presented study 
of the finite-temperature behaviour of ultra-cold Bose 
atoms in three-dimensional (3D) optical lattices by the 



slave fermion and the slave boson approaches to the Bose- 
Hubbard (BH) model. The finite-temperature phase dia- 
gram was also investigated by Gerbier^ in the context of 
ultra-cold bosons confined in optical lattice in the pres- 
ence of an additional potential. Three regimes can be 
recognised from the phase diagram: zero-temperature 
quantum phase, intermediate, where Mott insulator (MI) 
features persist but superfluid (SF) region is absent and 
the thermal region, where the Mott insulator properties 
disappear. Despite of the few theoretical approaches to 
the problem of the strongly interacting bosons at finite- 
temperatures many questions still remain open and un- 
solved. Especially studies of the phase transitions in the 
strongly correlated regime are scarce, where the repulsive 
energy is the main energy scale in the system and we are 
far from the limit of weakly interacting bosons. For these 
reasons, there is still a strong need for approximate but 
robust treatments of strongly correlated bosonic models 
in order to include the finite-temperature properties es- 
pecially relevant for phase diagrams for the 3D system 
where temperature-induced phase transition exists. The 
purpose of this paper is to present a robust theoretical de- 
scription of correlated bosonic systems which fulfils these 
goals. Our main idea is to focus on the degrees of freedom 
associated to the relevant physical variable associated to 
the Mott-superfluid transition, namely a quantum phase 
field U(l) rotor field, which is dual to the local occupa- 
tion number, and acquires dynamic significance from the 
boson-boson interaction^. 



II. THE MODEL 

In experimental parameter regime the bosonic atoms 
with repulsive interactions in a periodic lattice potential 
are perfectly described by a Bose-Hubbard model which 
is the simplest nontrivial model describing a bosonic 
many body system on a lattice which can not be mapped 
onto a single particle problem. Nevertheless it captures 
essential effects like a quantum phase transition from a 
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superfluid state to a Mott insulating state. The Hamil- 
tonian of the model reads 
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where, a] and stand for the bosonic creation and anni- 
hilation operators that obey the canonical commutation 
relations [a^aj] = 5ij, where rii = a\ai is the boson 
number operator on the site i. Here, (i,j) identifies sum- 
mation over the nearest-neighbour sites. Furthermore, 
tij is the hopping matrix element and describes the tun- 
nelling of bosons between neighbouring potential wells 
in the simple cubic lattice and Jjl/U = ///[/ + 1/2 is 
the shifted reduced chemical potential which controls the 
number of bosons, U > is the on-site repulsion. Due 
to the short range of the interactions compared to the 
lattice spacing, the interaction energy is well described 
by this term, which characterises a purely on-site inter- 
action. The interaction term tends to localise atoms to 
lattice sites. When the potential depth of the optical 
lattice is increased, the tunnelling barrier between neigh- 
bouring lattice sites is raised and the tunnelling matrix 
element t decreases. The on-site interaction U on the 
other hand is increased due to a tighter confinement of 
the wave function for bosons on a lattice site. Therefore 
the ratio U/t can be continuously adjusted over a wide 
range by changing the strength of the lattice potential. 
Finally we comment on the validity of the Hamiltonian in 
Eq. ([I]) . The Bose-Hubbard model can be obtained from 
many-body Hamiltonian with pseudo-potential interac- 
tion. However, we must assume that thermal and mean 
interaction energy are much smaller than the separation 
to the first excited band. Therefore, the exact value of 
the temperature (in the limit where the recoil energy is 
much smaller than maximum value of the lattice depth 
Er <C Vb) and can be calculated from the expression 
ksT <C Hujq where Hu)q is the energy separated a num- 
ber of vibrational levels in an infinite periodic potential. 
The energy to comparison can be taken from single site 
of the lattice. Moreover the hopping matrix elements are 
nonzero or non- negligible only to nearest-neighbours, so 
there are always negative in the lowest band. 



III. PHASE ACTION AND ORDER 
PARAMETER 



The functional integral representation of models for 
correlated bosons allows us to implement efficiently the 
method of treatment. The partition function is written 
in the form 
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and the bosonic path integral is taken over the complex 
fields di (r) with the action S given by 
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where /3 = l/k&T and T is the temperature. Since 
Hamiltonian is not quadratic in the fields a$ we have to 
decouple first the interaction term in Eq. by means of 
a Gaussian integration over the auxiliary scalar potential 
fields Vi (r) which periodic part V? (r) couples to the lo- 
cal particle number through the Josephson-like relation 
4>i ( r ) = ^ P ( r ) wnere fa ( T ) = dfa (r) I dr. The phase 
field satisfies the periodicity condition fa (/3) = fa (0) as 
a consequence of the periodic properties of the Vf (r) 
field. Next, we perform the local gauge transformation 
to the new bosonic variables 



Q>i (r) = bi (r) exp [ifa (r)] . 
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Using such a description is justified by the definition of 
the order parameter 
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(r)) = (bi (r) exp [ifc (r)]) = b ^ B (5) 



which non-vanishing value signals a macroscopic quan- 
tum phase coherence (in our case we identify it as SF 
state). The system can be then described by a macro- 
scopic wave function since the many-body state is a 
product over identical single particle states. Therefore 
a macroscopic phase is well defined on each lattice site 
and the system is superfluid. On the other hand the 
atom number per site is uncertain , and therefore one 
would find a random atom number in a measurement. 
In the large U limit the amplitude bo = (bi) (see Fig. 
[I]) has a nonzero value, but to achieve the superfluid- 
ity, the phase variables must also become stiff and, in 
consequence, i/jb = exp [^(r)] ^ 0- Furthermore, we 
parametrise the boson fields bi (r) = bo + b { (r) and re- 
strict our calculations to the phase fluctuations dropping 
the amplitude dependence^ which is justified in the large 
U/t limit. The U (1) group governing the phase field is 
compact and <\> (r) has the topology of a circle, so that in- 
stanton effects can arise due to non-homotopic mappings 
of the configuration space onto the gauge group U(l). 
Accordingly, the path integral reads 



{mi} 



f 4 

(0) / 



r),+27i 



d&(0)/ [D&(t)}... 

(6) 

which is performed by taking phase configurations that 
satisfy boundary condition fa (/?) — fa (0) = 2i\mi (m^ = 
0, ±1, ±2, . . .), where the winding numbers rii label the 
distinct homotopy classes of the U (1) group. Thus the 
paths can be divided into topologically distinct classes, 
characterised by a winding number defined as the net 
number of times the world line wraps around the system 
in the "imaginary time" direction. Integrating the action 
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Figure 1: The superfluid order parameter i/jb measuring the 
degree of the phase coherence as a function of the inter- 
action strength t/U for several values of the temperature 
k B T/U = 0.00,0.05,0.10,0.20 from the left to the right and 
fixed chemical potential. The dash-dotted line is the corre- 
sponding amplitude amplitude bo (see, Eq[5]). 



m Eq. over the bosonic fields we obtain the effective 
Lagrangian in terms of the phase-only variables 
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where the ampli- 



with the phase stiffnesses Ji • 

tude 6q = j) Uj + P^j /U originates from the saddle 
point condition 



dS[b, b] 



db 



(8) 



b=b 



The result of the gauge transformations is that we have 
managed to cast the strongly correlated problem into a 
system of mutually noninteracting effective bosons, sub- 
merged in the bath of strongly fluctuating U (1) phase 
fields, whose dynamics is governed by the energy scale set 
by the on-site interaction U that drives the Mott transi- 
tion. 



IV. TREATMENT OF THE ACTION OF 
QUANTUM ROTORS 

Now, we devise a systematic way of treatment for the 
fluctuating phase fields contained in the action in Eq. ^ 
that enables us to obtain an effective non-linear sigma- 
field theory that respects the symmetry properties of the 
model and satisfies the Mermin- Wagner theorem, thereby 
improving the pure mean-filed approach known from its 
restricted ability to deal with the spatial fluctuations. To 



proceed, it is convenient to replace the phase degrees of 
freedom by the complex field ipi = e^ i<Kj ^ which satisfies 
the periodic boundary condition ipi (/3) = ipi (0). This 
can be done by implemented the Fadeev-Popov method 
with the Dirac delta functional resolution of the unityP^ 
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where we take ipi as continuous variable but constrained 
(on the average) to have the unimodular value. We can 
solve the constraint by introducing the Lagrange multi- 
plier A which adds the quadratic terms (in the ipi fields) 
to the action Eq. The partition function is written 
in form 
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where the free energy per site T — — In Z/f3N is given 
by: 
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and jij (r, r') = (exp{— i (r) — <pj (r')]}) is the two- 
point phase correlator associated with the order parame- 
ter field, where (...) is the averaging with respect to the 
action in Eq. 

The action with the the topological contribution, after 
Fourier transform, we write as 



(12) 



where T^ 1 (uji) = X — Jk + 7 -1 (uje) is the inverse of the 
propagator and uj£ = 2tt£/P (£ — 0, ±1, ±2, . . .) stands 
for the Bose-Matsubara frequency, while the phase cor- 
relator, after Fourier transform, can be written as: 
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tition function for the set of non-interacting quantum 
rotors. Within the phase coherent state the order pa- 
rameter is given by 
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Figure 2: Phase diagram for BH model on a cubic lattice as 
a function of chemical potential and temperature. The Mott 
insulator phase is found only at T — within each lobe of inte- 
ger boson density since for T=0 this state is incompressible.; 
at T > a region with finite compressibility emerges (due 
to the thermal activation), so we have here the "disordered 
state" without phase coherence. Above the critical surface 
the superfluid region takes place. Top panel: the projection 
of the critical surface onto the plane temperature-chemical 
potential in a form of the density plot. 



For the simple cubic lattice we write Jk = (12t + p) t^/U 
with the dispersion = 2t (cos k\ + cos k^ + cos k^). 
The phase boundary is determined by the divergence 
of the order parameter susceptibility Tk=o (^£=o) = 00 ? 
which determines the critical value of the Lagrange pa- 
rameter A = Ao that stays constant in the whole ordered 
phase. After summation over Matsubara frequency the 
superfluid state order parameter becomes 



coth 



(15) 



In the above equation A£ = (J — Jk) /U + v 2 (/jl/U) and 
v (fi/U) = frac (fi/U) — 1/2, where frac (x) = x—[x] is the 
fractional part of the number and [x] is the floor function 
which gives the greatest integer less then or equal to x. 
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Figure 3: Bosonic occupation number asa function of the 
chemical potential, for several values of the temperature as 
indicated in the plot calculated for the extreme interaction 
limit (U/t = oo, to highlight the sole temperature effect). 



V. RESULTS 

When the strength of the interaction term relative to 
the tunnelling term in the Bose-Hubbard Hamiltonian 
is changed, the system reaches a quantum critical point 
in the ratio of U/t, for which the system will undergo 
a quantum phase transition from the superfluid ground 
state to the Mott insulator ground state. In three dimen- 
sions, according to Mermin- Wagner theorem this phase 
transition can occur also at non-zero temperature. The 
finite-temperature phase diagram of the model can be 
calculated from Eq. (15) by introducing the density of 



states for simple cubic lattice in order to perform the sum 
over the lattice wave vectors: 



p(0 = 




(16) 



with a\ = min (—1, —2 — £/t) and = max (1, 2 — £/£); 
K (x) is the elliptic function of the first kincP^. 

A lobe-like structure (see Fig. |2]), similar to the zero- 
temperature case, becomes flat with increasing temper- 
ature. The lobes with a larger boson occupation num- 
ber are more stable against temperature. The stabil- 
ity comes from higher values of the repulsive energy U. 
Therefore, at temperature T = the interaction in the 
system^* governs the quantum phase transition. By de- 
creasing the value of the repulsive energy we can achieve 
superfluid phase. In real physical realizations of the BH 
model thermal excitations are always present and also 
can activate a phase transition. We see it clearer by 
calculating a bosonic occupation number. Decompos- 
ing the phase field in terms of a periodic field and lin- 
ear in r term we calculate the effects of the fixed bo- 
son number ub = N~ x (di (r) (r)) in the system. 
The total boson density n# = + Sn^ consists of the 
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Table I: Comparison of the maximum of the critical value 
for t/U parameter (as a function of the normalised chemi- 
cal potential /jl/U) at the tip of the first (ns = 1) MI lobe 
for square lattice with several numerical (QMC - quantum 
Monte- Car lcP^, DPT - diagrammatic perturbation theor^P^) 
and analytical work: SCPT - strong coupling perturbation 
theory 20 . QRA - our calculations using quantum rotor ap- 
proach) . 



0.20 

k B T/U 





QMC 


DPT 


PA 


QRA 


t/U 


0.03408(2) 


0.03407 


0.034737 


0.03215 


fJL/U 


0.389 


0.393 


0.37905 


0.41 



occupation number for neutral bosons and a contri- 
bution 5ni) from a fluctuating phase field. For T = 
we recognise a steps of fixed integer filling of bosons (see 
Fig. |3j. With increasing temperature, typical for the 
Mott state, steps-like profile becomes smoother. There- 
fore, bosons placed in the Mott state get energy required 
to move from one lattice site to another from thermal 
fluctuations. The temperature k^T /U ~ 0.2, where the 
occupation number characteristic becomes flat, is simi- 
lar to recognised as a melting temperature for the con- 
densate slowly loaded into the optical potential in the 
presence of an smoothly varying trap^. The increas- 
ing value of the thermal energy, in analogy to the de- 
creasing repulsive interaction^, can lead to the situation 
where the sharp steps of the MI state become indistinct 
(see Fig. [3]). The Mott-insulator to superfluid quan- 
tum phase transition is rigorously present only at zero 
temperature, whereas at finite-temperature thermal fluc- 
tuations induce a phase transition between a superfluid 
and a normal phase. However, at sufficiently low tem- 
peratures, a remnant of the insulating phase still persists 
within the normal phase. In these conditions it is possi- 
ble to observe a sharp crossover between a compressible 
normal fluid and a phase characterised by a vanishing 
compressibility (see, Fig. [4|. Regarding the comparison 
of our method with the previous approaches we found 
that our results are in good agreement with other calcu- 
lations including numerical quantum Monte-CarlcP^, di- 
agrammatic perturbation theory 19 ) and analytical works 
based on the strong coupling perturbation theory^, see 
Table [T] and Fig. [5] 



VI. FINAL REMARKS AND OUTLOOK 

To conclude, in this paper we have presented a study of 
the finite-temperature transition of the three-dimensional 
Bose-Hubbard model relevant for Bose-Einstein conden- 
sates trapped in three-dimensional optical lattice poten- 
tials which allow to enter a new regime in the many body 
physics of ultracold atomic gases. Ultracold atoms in op- 
tical lattice potentials represent a rather perfect realiza- 
tion of the Bose-Hubbard model with a high degree of 
control. Parameters like the ratio between the on-site in- 
teraction and tunnel coupling or the filling factor can be 
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Figure 4: Density plot of the normalised compressibility k — 
k/U, where k = dns/d/n calculated numerically as a function 
of temperature and chemical potential for t/U = 0.012. Dark 
(red) shading around the integer values of ji/U corresponds 
to high values of k whereas the light (blue-to-green) shading 
marks the region of diminishing compressibility. 
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Figure 5: The comparison of the zero-temperature phase dia- 
gram calculated from the diagrammatic perturbation theor^P^ 
(DPT) and our results obtained in the frame of quantum ro- 
tor approach (QRA) with ub = 1, 2, 3 and 10. We found 
the critical value (t/U) CTit (see also Table fi| for the tip of the 
nth Mott-insulator lobe is always slightly lower than obtained 
from DPT and Monte-Carlo simulations^. 



widely adjusted, and observables like long range phase co- 
herence or number statistics can be measured. In order 
to approach these issues theoretically we employed the 
U(l) quantum rotor approach and a path integral for- 
mulation of quantum mechanics including a summation 
over a topological charge, explicitly tailored for the BH 
Hamiltonian. This method can give the thermodynamics 
of the Bose-Hubbard model in the limit of strong interac- 
tions. Our aim was then to analyse the phase transitions 
that may occur in such system at finite-temperature, and 
determine the general features of the associated phase di- 
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agrams. We demonstrated the evolution of zero temper- 
ature Mott lobes and Mott plateaus when the tempera- 
ture is increasing. The technique used in this paper can 
be conveniently extended to more general situations, in- 
cluding i.e. multi-species bosonic systems. Hamiltonians 
other than the pure Bose-Hubbard Hamiltonian could be 
realized and studied by the method outlined in three- 
dimensional optical lattice potentials. For example by 
using a multi component gas and interspecies Feshbach 
resonances, an intriguing system could be created. Fur- 
thermore, it should be possible to consider other inter- 



actions as effective next neighbour interaction, dipolar 
interactions or spin interaction due to spin dependent 
tunnelling. Other generalisations of the on site Hubbard 
model are of course possible by including the influence of 
the disorder in the description. 
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